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We study the classical non-equilibrium statistical mechanics of scalar field theory on the lat- 
tice. Steady states are analyzed near and far from equilibrium. The bulk thermal conductivity is 
computed, including its temperature dependence. We examine the validity of linear response pre- 
dictions, as well as properties of the non-equilibrium steady state. We find that the linear response 
theory applies to visibly curved temperature profiles as long as the thermal gradients are not too 
strong. We also examine the transition from local equilibrium to local non-equilibrium. 
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The understanding of the dynamics of non-equilibrium field theories is important to many areas in physical sciences, 
from processes in inflation or baryogenesis in the early universe, transport processes in condensed matter, to the 
possible states of hadronic matter in heavy ion collisions, such as quark-gluon plasma, disoriented chiral condensates 
and color superconducting states. However, many problems linger even in the basic understanding of non-equilibrium 
statistical mechanics and transport 0| . It is clearly desirable to address these problems regarding the non-equilibrium 
dynamics of field theories while making as few assumptions about the dynamics of the theory as possible. In this 
work, we study the steady state dynamics of classical massless (j>* lattice field theory in {d + 1) dimensions {d = 1, 3), 
under weak and strong thermal gradients. We study the classical field theory on a lattice since the problems are 
well defined and techniques are available to construct non-equilibrium states from first principles. It is not clear to 
us how to address the questions we pose in the full quantum theory without making some drastic assumptions. Our 
computations are non-perturbative. Furthermore, not only is the temperature within the system dynamical, but 
even the question of whether local equilibrium is achieved is not an assumption but is determined dynamically by 
the system. The equilibrium properties of classical 0^ theory have been studied in the past, including its ergodic 
properties 1^,^, Hamiltonian dynamics and phase transitions Q|. However, the kind of questions we address here have 
not been answered in the previous literature. 

Near equilibrium, linear response theory supposedly holds; yet there is no means to address its regime of validity 
within the theory itself since the computation is performed in equilibrium. By explicitly constructing non-equilibrium 
steady states near equilibrium, we examine the validity of the linear response theory. We then construct steady states of 
the system under stronger thermal gradients and study their physical properties. Here, fundamental questions arise, 
such as under which conditions local equilibrium is achieved. Local equilibrium, an assumption that equilibrium 
concepts can be applied locally to a problem which might be globally non-equilibrium, is widely used In fact, 
thermalization is a concept which is assumed, often tacitly, in many applications of non-equilibrium physics. However, 
as a system moves away from equilibrium, what precisely constitutes 'local equilibrium' becomes unclear and it is of 
interest to understand what kind of deviations develop from it and why. These problems are also non-trivial from the 
statistical mechanical point of view; even small departures from equilibrium into non-equilibrium steady states (such 
as the those we study) already lead to peculiar behaviors including a divergent Gibbs entropy Sg (Sg — oo) and a 
multi- fractal steady state measure [^,|l| . The steady-state distributions far from equilibrium are not well understood 
classically and very little is known concerning their quantum counterparts. We would like to see to what extent these 
peculiar statistical measures influence the steady state thermal profiles, T{x). 

Classical field theories are relevant to the high temperature dynamics of quantum field theories and have proved 
effective, for instance, in computing finite temperature properties of the standard model |0,||j8|. Furthermore, non- 
equilibrium dynamics of classical field theories is of interest in its own right, an understanding of which is essential to 
understanding the dynamics of the quantum theory. The approach we adopt here can be applied to classical lattice 
theories quite generally. Transport properties have been previously studied in scalar quantum field theory using linear 
response theory |^ and the theory is known to have a classical, finite temperature limit for correlation functions [^ . 
Yet the question of how to relate our results to those results is far from trivial and will not be pursued here. 

We start with the Lagrangian 

1 / f)M;^\\ 7.2 

(1) 

This model, when discretized, reduces to a model of lattice vibrations with quartic anharmonicity, with the following 
dimensionless Hamiltonian 




+ (V0r)' + Ut 



(2) 



Here tt ~ d(f)/dt, r runs over all sites in the lattice, and the lattice derivative has components ^k4'r = '/'r+ek ~ 4'r 
(ck is the unit lattice vector in the fc-th direction). The two theories are related by discretization and the rescalings, 
4>r{t) = ag(f){r,t), t = t/a, r = f/a, where a is the lattice spacing. The equations of motion, □(/> — —(j)^, are solved 
on a spatial grid, using two methods: fifth and sixth order Runge-Kutta, and leap-frog algorithms |jll|. 

In order to generate a stationary non-equilibrium statistical ensemble, thermal boundary conditions are imposed 
on the equations of motion. Specifically, at a; = and x = L, we add two time-reversal invariant fields which 
act to dynamically thermalize these boundaries at given temperatures Ti and T2 [jl^ (a more detailed account will 
be given elsewhere). For d > 1 we impose periodic boundary conditions on the other directions. Apart from the 
thermal boundary conditions, the system evolves according to the dynamics dictated by the Hamiltonian (^). We 
used from 10^ to 10^ time steps of dt from 0.1 to 0.001, with observables being sampled every Ai = 20 ~ 100 di. 
In d = 1, the lattice size was varied from L = 20 to 8000, while in d = 3 it ranged from 50 x iV x iV {N ranging 
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from 3-20) to 1000 x 3 x 3. We have verified that when Ti = T2, these boundary conditions dynamically set all the 
temperatures inside the system to be equal to the boundary temperatures and reproduce the equilibrium canonical 
measure Peq{'^T4') exp[— 77(7r, (^)/ri] at all points. 

By controlling Ti and T2 we can begin to explore the non-equilibrium steady state. One question we would like to 
address is how the temperature profile T[x) behaves. Near equilibrium one would expect a linear profile, but beyond 
that the shape is unknown. In our near equilibrium simulations, Ti < T2, we find a linear temperature profile and 
recover transport given by Fourier's law, as shown in Fig. 1 (a). However, far from thermal equilibrium (Ti <C 12), 
the temperature profile develops significant curvature, seen in 1 (b),(c) for the ratio T2/T1 — 10, 20. It is of interest 
to understand the physics behind these temperature profiles. We would further like to understand until what point 
linear response and local equilibrium provide reasonable descriptions. 

Equilibrium: Ti = T2. A standard approach to thermal conductivity utilizes equilibrium correlation functions to 
compute near-equilibrium transport. This linear response approach uses the Green-Kubo formula, 

k{T) = Jdr (rO-(r,i)r°-(ro,0))^^, (3) 

where the autocorrelation function is evaluated in the canonical ensemble, Ti — T2. For our lattice calculation, we 
replace / dr with a lattice sum. It is interesting to note that the integrand in (3) has been argued to develop a long 
time tail behavior of ^ t^'^/'^, leading to the divergence of (3) in d = 1 |l^. In Fig. 2 (top), a typical autocorrelation 
function for (H-l) dimensions is plotted to several hundred times the mean free time, which is well into the regime 
where long-time tails would be evident. The time integral is given in Fig. 2 (bottom), showing that the integral (3) 
is finite, which we attribute to the 'on-site' nature of the (j)^ interaction, in contrast to some of the other models p4[ . 
Interestingly enough, we do find that the transient behavior of the Green-Kubo integrand is quite close to t^^/"^ up 
to a few ten times the mean free time, after which it decays much faster. Consistent results are found in d = 3 as 
well. The computed n{T) is shown in Fig. 3 (top). 

JVear Equilibrium: Ti < T'2. Near equilibrium, we find that a linear temperature profile emerges dynamically. The 
thermal conductivity k is then obtained through Fourier's law: 

niT) = = "a*0V0, (4) 

where {T'^^)„e is the heat flux averaged over the non-equilibrium steady state. Attention is paid to verifying the 
linear response properties by varying the temperature difference \T2 — Ti \ around the same average temperature. T{x) 
is the local temperature defined through an ideal gas thermometer, by T(x) = {tt'^ (x)) ^ e , where Tr(x) is the momen- 
tum density. This will serve as a convenient definition as long as local equilibrium is achieved and the momentum 
distributions are gaussian. Here (■ • indicates the ensemble average over the non-equilibrium steady state. To 
obtain the transport properties, each simulation is run long enough for observables such as T'^^, the energy density 
as well as distribution functions to converge. 

In Fig. 3 (top), we compile the Green-Kubo and direct measurements of k, plotted as a function of T. We find that 
these independent computations are quite consistent with each other, k is found to have a temperature dependence 

_ Arp-y 17 = 1.38(2), = 2.72(4) (1-1-1) dimensions 
^ ' 1 7 = 1.64(4), A = 9.1(2) (3+1) dimensions ' 

This behavior is similar to that of lattice phonons at high temperature ||l^] . We have also verified that a sensible bulk 
behavior exists, as shown in Fig. 3 (bottom); the thermal conductivity is independent of L when it is larger than the 
mean free path, which, on the lattice, is of order of the conductivity. 

In trying to understand near equilibrium physics, one might be tempted to assign a statistical measure, such as 
Pne{t^, <P) ~ exp[— i7(7r, 0)/T(x)] to the non-equilibrium stationary state, or similar measures which assume some form 
for T{x). Strictly speaking, this is not correct; the phase space measures which describe steady state non-equilibrium 
systems are multi- fractal (whether one studies shearing, heat fiow and so forth), converging to smooth distributions 
(Boltzmann) only in the equilibrium limit |^jl] . An important consequence is that the dynamical space is necessarily 
of lesser dimension than the equilibrium phase space. This in turn results in additional correlations of 7r(a;) and 4){x). 
One can also see that the non-equilibrium measure is not locally Boltzmann since quantities such as {'K{x)(t>{x')) / 
for X ^ x' which is also reflected in the non-zero energy flow, while (7r(x)) = near and far from equilibrium. 

Far From Equilibrium: Ti ^ T2. When the temperature gradients are larger and we are no longer in the linear 
regime, the temperature profile becomes visibly curved. An example of such a non-linear temperature profile are 
given in Fig. 1. It should be noted that inside the boundaries, the dynamics is that of only the 0"* theory and the 
temperature profile is determined by it. When the temperature varies substantially in the system, one cause for the 
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non-linearity is the temperature dependence of the thermal conductivity. If this were the only cause of non-linearity, 
the temperature profile can be determined within the region by integrating Eq. (4) (when 7 7^ 1) 



T{x) = Ti 



(6) 



T{x) is a function only of x/L so that it is consistent with a smooth contimmm limit. This expression for T{x) 
provides very good descriptions of the measured profiles with significant curvature, as can be seen in Fig. 1 where 
Eq. (6) (dashes) is almost indistinguishable from the measured steady state profiles (solid) . This indicates that linear 
response extends well beyond the regimes of small temperature differences, when applied locally. A priori, this was 
not clear. 

By continuing to increase the ratio T2IT1 in the simulation, we do eventually reach steady state situations where 
the linear response formula no longer works. At this point, we begin to see indications that the concept of local 
equilibrium also becomes more tenuous. To analyze these questions more concretely, we need a dimensionless measure 
of how strong the thermal gradient is. A natural choice we adopt is XVT/T, where A is the mean free path. In our 
model, the heat capacity per unit volume, Cy, and the sound speed, c^, are of order unity, so that elementary kinetic 
theory suggests that the mean free path is A ~ dn, where d is the spatial dimension. To quantify the departures from 
the linear response formula, we plot in Fig. 4 (for one spatial dimension) the deviation of the measured heat flux to 
that obtained from linear response prediction using Eqs. (^) — (^ (denoted {T°^)ne and {T'^^)Lji, respectively), as a 
function of the quantity k(T)WT/T. The figure includes different lattice sizes and temperatures. We see that for 



VT 
<T)^ « 1 



(7) 



linear response theory holds quite well. This includes systems with significant curvature in T{x) such as Fig. 1. Even- 
tually, for sufficiently strong gradients, the measured heat flow begins to deviate from the linear response results; the 
system does not conduct heat as well as its linear response theory prediction. At this point, we observe simultaneously 
the departure of other quantities from a 'local equilibrium' characterized by gaussian momentum distributions. To 
see the departure from local equilibrium, we follow the behavior of various observables, which include the momen- 
tum cumulants, the steady-state momentum distributions /(tt^) as well as heat flux and correlation functions. For 
instance, in equilibrium, (n^ (x)) / {tt"^ (x))'^ = 3, {n^ (x)) / {n^ (x))^ = 15 and so forth. In the regime where the Hnear 
response theory breaks down, the momentum distributions become more sharply peaked and are no longer gaussian 
even in the steady state, and the deviations in the cumulants become apparent. This indicates that at least in our 
theory, higher order corrections to linear response are not that well founded since the concept of temperature becomes 
tenuous at this point. Since the heat flow is a constant, we also note that the local equilibrium condition (|^) is more 
likely to be satisfied at the higher temperature end. Similar behavior is seen in (3-1-1) dimensions. 

We have further examined the behavior of the (coarse grained) Boltzmann entropy Sb- While Gibbs entropy is 
known to be singular, any similar divergence in the coarse grained Sb would only be evident if one extrapolated 
the measured values to the continuum limit. Instead, we consider if the notion of local entropy changes significantly 
as the system moves away from equilibrium. To this end, we compute the Boltzmann entropy Sb from the 1- 
and 2-body densities f'^^\Tr{x),(j){x)) and /'^^^(7r(a;), ^(x), 7r(x'), (/)(a;')) in the non-equilibrium steady states, from 

— — J dfi^''^ f^^Hog/'^'^K We find that S^^^ does not shift noticeably from its equifibrium value regardless of 
how far the system is from equilibrium. Further, 5^^^ (< 25''^^^) is only slightly less than its upper limit 25^^^ and 
remains so even far from equilibrium. So unlike Sg ^< ^^b^^^ is found to be rather insensitive to the non- 
equilibrium nature of the system. While this could be a manifestation of coarse graining, it does suggest that some 
local thermodynamic concepts might still be useful in making connections with non-equilibrium thermodynamics. 

The behavior of scalar lattice field theory near and far from equilibrium has been explored, and the present study 
allows us to establish a number of features regarding the non-equilibrium stationary state. The bulk thermal conduc- 
tivity and its dependence on the temperature over few decades was found. To our knowledge, such a computation 
from first principles, without assuming linear response, has not been performed previously. This was done using both 
linear response (equilibrium) and direct (near equilibrium) approaches. The d = 1 Green-Kubo integrals in our theory 
are non-divergent and agree with the direct computation, in contrast to some of the other models. We find that linear 
response theory is quite robust and works well even when the steady state thermal profile has significant curvature, 
and as a consequence we derive an analytic description for the profile. By driving the system farther from equilibrium, 
we are able to see when linear response breaks down. Surprisingly, this is found to be near the same point where 
local equilibrium becomes noticeably violated. A sufficient condition for noticeable departures from linear response 
and local equilibrium in d = 1,3 seems to be AVT/T > 1/10. It would be interesting to examine the dynamics of 
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non-equilibrium phase transitions or explore dynamical cooling of boundary temperatures to provide a means to access 
more complex non-equilibrium environments, and in particular, the dynamics of ultrarelativistic heavy-ion collisions. 
We acknowledge support through grants at Keio University and DOE grant DE-FG02-91ER40608. 
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FIG. 1. (a) Linear thermal profile (solid) near equilibrium for L=800, d = 1 compared to a linear fit (dashes). The 
temperature was sampled at At = 1 for 5 x 10® points, (b) Sagging profile for a non-equilibrium steady state in (1+1) 
dimensions far from equilibrium (solid) compared to theoretical predictions (dashes), where (Ti, T2) = (0.1, 1). The temperature 
was sampled 10® times every At = 0.25. (c) A curved profile in (3+1) dimensions (solid) with the theoretical fit (dashes) for 
a 1000 X 3 X 3 system with (TijTg) = (0.5, 10). The temperature was sampled 4 x 10® times every At = 0.1. In both (b) and 
(c), the fits are virtually indistinguishable from the temperature profiles. 
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FIG. 2. (Top) Time dependence of \C{t)\, where C{t) = f da (T°"^(a;, t)T°"^(a;o, 0))^^ /T^, for L = 100, up to a time 

t — 2x 10*, well into the region where long-time tails should be evident. Inset: Short time decorrelation. (Bottom) Green-Kubo 
integral up to time t, converging to k in 1+1 dimensions. 
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FIG. 3. (Top) Thermal conductivity k obtained from direct (□ 1 + 1-d; o 3+1-d) and Green-Kubo (x) measurements for var- 
ious lattice sizes L, and the power law fit (dashes, dots). (Bottom) L dependence of k indicating bulk behavior for temperatures 
(upper to lower) T = 1/10, 1/4, 1 (solid, l+l-d) and T = 1 (dashes, 3+1-d). 



2 




0.0 0.1 0.2 0.3 



/c(T)VT/T 

FIG. 4. Breakdown of the linear response theory: The deviation of the measured energy flow {T^'')ne from its linear 
response expectation, {T'^^)lr, is shown as a function of kVT/T for d = 1 non-equilibrium steady states. As discussed in the 
text, local equilibrium fails simultaneously. 
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